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1. Introduction 

The identification of modular properties in integrable tlieories lias always represented a 
progress in our understanding of the underlying physical and mathematical structures of 
the model. Modular transformations can appear in a physical model in many different 
ways. For example, in two-dimensional Conformal Field Theory (CFT) the modular 
invariance of partition functions on a toroidal geometry [1] has lead to the classification 
of entire classes of theories [2] and, in many cases, has opened the way to their higher 
genus investigation, which is a very important tool in applications ranging from strings 
to condensed matter physics. In a different context, modular transformations are at the 
root of electromagnetic duality in (supersymmetric) gauge theories |3|. In this case the 
modular transformations act on the space of parameters of the theory, rather than as 
space-time symmetries. 

In integrable two-dimensional lattice models modular properties were crucial in 
a Corner Transfer Matrix (CTM) approach [4J to get the so called highest weight 
probabilities |5l [6]. More recently, the problem of computing bipartite entanglement 
Renyi entropies (or their Von Neumann limit) in integrable spin chains has been related 
to the diagonalization of the CTM of the corresponding 2D lattice integrable model. 
The expressions for these quantities often present a form that can be written in terms 
of modular functions, like the Jacobi theta functions, suggesting links with a modular 
theory in the space of parameters. In the XY chain, for example, modular expressions 
have been found for the entanglement entropy in [7]. Similar modular expressions have 
been obtained in the calculation |8], [9] of entanglement entropies in the XYZ spin chain, 
notoriously linked in its integrability to the 8- vertex lattice model The XYZ model 
is one of the most widely studied examples of a fully-interacting but still integrable 
theory, playing the role of a sort of "father" of many other very interesting and useful 
spin chain models, like the Heisenberg model and its XXZ generalization. Here we 
would like to start the investigation of its modular symmetries, which - let us stress 
it again - are not spacetime symmetries, like in CFT, but symmetries in the space of 
parameters. We believe that the exploiting of modular properties hidden behind the 
Baxter reparametrization of the XYZ chain can lead, similarly to the aforementioned 
examples of conformal and gauge theories, to new tools in the investigation of crucial 
properties of the XYZ chain and its family of related models. 

In this paper we show that the XYZ Hamiltonian is invariant under a set of 
transformations that implement a modular group covering the full plane of parameters 
of the model in a very specific way. This, in our opinion, prepares the way to the 
application of the modular group to the study of more complicated quantities, that 
we plan to investigate in the near future. For instance, this modular invariance may 
play an important role in a deeper understanding of the Renyi and Von Neumann 
bipartite entanglement entropies, which are well known to be excellent indicators to 
detect quantum phase transitions. The paper is organized as follows. In section 2 
we recall basic facts about the Hamiltonian formulation of the XYZ spin chain, the 
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properties of its phase space diagram and the symmetry transformations leading one 
region of the diagram into another. In section 3 we recall the Baxter solution allowing 
to relate different portions of the phase space by transformations of the parametrization 
of the coupling constants of the model in terms of elliptic parameters. Section 4 is 
devoted to the analytical extension of Baxter solution by switching to new parameters 
that make the modular properties more evident. This shows how one of the phase 
space regions can be mapped into the others by applying a modular transformation. 
The generators of these modular transformations are identified in the next section 5, 
where the characterization of the actions of the full modular group is described. Our 
conclusions and outlooks are drawn in section 6. 



2. The XYZ chain and its symmetries 

The quantum spin-| ferromagnetic XYZ chain can be described by the following 
Hamiltonian 

n 

where the a" (a = x, z) are the Pauli matrices acting on the site n, the sum ranges 
over all sites n of the chain and the constants J^-, Jy and take into account the degree 
of anisotropy of the model. We assume periodic boundary conditions. 

In figurejljwe draw a cartoon of the phase diagram of the XYZ model in the ^ j^, 'j- 
plane. We divide it into 12 regions, whose role will become clear as we proceed, named 
Ia,b,c,; IIa,b,c,; Hla.b.c, • BoM (red-onlinc) continuous lines indicate the gapless phases. To 
understand better the model, let us first look at the Jy = Jx line. Here we have the 



familiar XXZ model and we recognize the critical (paramagnetic) phase for 



< 1 and 



the ( ant i-) ferromagnetic Ising phases for 



> 1. The same physics can be observed 



for Jy = — Jx- There, one can rotate every other spin by 180 degrees around the x-axis 
to recover a traditional XXZ model. Note, however, that J^ changes sign under this 
transformation and therefore the ferromagnetic and anti-ferromagnetic Ising phases are 
reversed. The lines Jz = ±Jx correspond to a XYX model, i.e. a rotated XXZ model. 
Thus the phases are the same as before. Finally, along the diagonals Jy = ±Jz we have 
a XYY model of the form 

HxYY = [J^«+i + Jy (o-^o-^+i ± «+i)] • (2) 

n 

Thus the paramagnetic phase is for ^ > 1 and the Ising phases for ^ < 1, with 
the plus sign {Jy = J^) for Ising ferromagnet and the minus {Jy = —Jz) for the anti- 
ferromagnet. 

In this phase diagram we find four tri-critical points at ^j^, = (±1,±1). As 
each paramagnetic phase can be described as a sine-Gordon theory with = Svr at the 
anti-ferromagnetic isotropic point, flowing to a = theory at the ferromagnetic 
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Figure 1. Phase Diagram of the XYZ model in the j^j plane. 

Heisenberg point, the critical nature of these tri-critical points is quite different. 
Assuming > 0, at Ci = ( J, = (-1,1) and C2 = ( J, = (1,-1) we 
have two conformal points dividing three equivalent (rotated) sine-Gordon theories with 
(3^ = 8n. At El = (^j^, = (1, 1) and E2 = J, = (-1, -1), we have two (3'^ = 
points which are no longer conformal, since the low energy excitations have a quadratic 
spectrum there. The former points correspond to an antiferromagnetic Heisenberg chain 
at the BKT transition, whilst the latter correspond to an Heisenberg ferromagnet at 
its first order phase transition. The different nature of the tri-critical points has been 
highlighted from an entanglement entropy point of view in [9j . 

We have seen that the model is invariant under a | rotation Rq, = Yin e'^ of every 
spin around one of the axes a = x,y, z, followed by the interchange of the couplings R/j,^ 
on the perpendicular plane (for instance, for a = x,we need Ry^^ Jy = Jz, i^y,z Jz = Jy)- 
Note that the composition of two rotations, say Rj, ■ R2 is compensated by the exchange 
of the coupling in opposite order: R^. y-R^, j.. Similarly, the inversion of every other spin 
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in a plane Pq, = Yl^ ctj^ can be compensated by changing the signs of the couphngs on 
that plane P^^ (where, P^ ^ Jy = —Jy, Py^z Jz = ^Jz)- Three operations are sufficient 
to generate all these symmetries, for instance R^., and P^; (and Ry,2, Ra;,^ and Py^z)- 
These generators relate the different regions of the phase diagram: R^Ia = Ifo, 
Ryla = Ilia, Pxia = Ic and so on. Note that the action of the Rq,,^ is to exchange the 
coupling, and we have six different orderings of the three coupling constants. The action 
of Py^z additionally doubles the possible phases given by a set of three couplings, giving 
the total of twelve regions of figure [1} We will show that these operators are equivalent 
to an extension of the modular group, acting in parameter space. 



3. Baxter's solution 

The Hamiltonian ([T| commutes with the transfer matrices of the zero-field eight vertex 
model |1] and this means that they can be diagonalized together and they share the 
same eigenvectors. Indeed, as shown by Sutherland [10], when the coupling constants 
of the XYZ model are related to the parameters F and A of the eight vertex mode||] at 
zero external field by the relations 

J, : J, : J, = 1 : r : A , (3) 

the row-to-row transfer matrix T of the latter model commutes with the Hamiltonian 
H of the former. In the principal regime of the eight vertex model it is customary |3] 
to parametrize the constants F and A in terms of elliptic functions 

1 + sn^(iA; A;) cn(iA; A;) dn(iA; /c) , , 

" 1 - A;sn2(iA;fc) ' " l-fcsn2(iA;fc) ' ^' 

where sn{z\ k), cn{z; k) and dn{z; k) are Jacoby elliptic functions of parameter k, while 
A and k are the argument and the parameter (respectively) whose natural regimes are 

0<A;<1, 0<A< K{k') , (5) 



K{k') being the complete elliptic integral of the ffist kind (A. 7) of argument k' = 

In light of (|3|, without loss of generality, we rescale the energy and set = 1. We 
recall that the parametrization ^ of the parameters is particularly suitable to describe 
an anti-ferroelectric phase of the eight vertex model, corresponding to A < —1, |r| < 1. 
However, using the symmetries of the model and the freedom under the rearrangement 
of parameters, it can be used in all cases by re-defining the relations that hold between 
A and F and the Boltzmann weights (for more details see [4]). We can apply the same 
procedure to extract the physical parameters of the XYZ model in the whole of the 

§ We adopt the conventions of [4 on the relation among F, A and the Boltzmann weights of the XYZ 
chain. 
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phase diagram in a fairly compact way. In fact, when \Jy\ < 1 and \Jz\ < 1 we have: 

IJ.. _ j^.l _ IJ.. _)_ j^.l 

r -- 

A -- 



1 Jz Jy 






\Jz Jy 


+ 


t/j -|- t/jyl 



(6) 



\Jz Jy \ ~^ 

while for I Jyl > 1 or I J^l > 1, we get: 



\Jz H" 



min 


1, 


\J.-Jy\-\Jz + Jy\ 




2 




max 


1, 


\Jz-Jv\-\Jz+Jy\ 




2 





4- ■ sgn IJ, - J„ 



\Jz + J 



I/I / ' 



(7) 



A 



1 

2. 



\Jz 






max 


1, 


IJz-JyI-IJ.+Jal 




2 





where F and A are given by (|4]). For instance, for \Jy\ < 1 and < —1 we recover 
Jy = r and Jz = A, while for Jy > 1 and jJ^I < 1 we have Jy = —A and Jz = — F. 
In figure [T] we show the phase diagram divided in the different regions where a given 
parametrization applies. 

While Baxter's procedure allows to access the full phase-diagram of the model, it 
introduces artificial discontinuities at the boundaries between the regions, as indicated 
by the max and min functions and the associated absolute values. In the next section 
we will show that, with a suitable analytical continuation of the elliptic parameters, 
we can extend the parametrization of a given region to the whole phase diagram in a 
continuous way, which numerically coincides with Baxter's prescription. 



4. Analytical extension of Baxter's solution 

To streamline the computation, it is convenient to perform a Landen transformation on 
Q and to switch to new parameters (u, I) 

k)X , 



/ = 



u = (l 



so 



thaf 



1 + ' 

^ _ 1 + A; sn2(iA; k) _ 1 

1 — k sn^ {iX; k) dn(iM; /) ' 
^ cn(iA; k) dn(iA; k) cn(iM; I) 
1 — k sn^liX] k) dn(iti; /) 

The natural domain of (A, k) corresponds to 

< M < 2K{1') < / < 1 . 

We also have the following identities 

1-k 



I' = v/f^T^ 



1 + A; ' 

See, for instance, table 8.152 of I 



K{1) = {l + k)K{k) 



Kil') 



1 + k 



(8) 

(9) 

(10) 

(11) 

K{k'), (12) 
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from which it follows that the Landen transformation doubles the elliptic parameter 



Tik) 



■K{k') 
' K(kj_ 



2i 



K(l) 



2r(/). 



Using (A.14)) we have 
/ = 



r2 



A2 _ r2 ■ 

Furthermore, elliptic functions can be inverted in elliptic integrals and from (|10 
have 



(13) 



we 



m 



^{l-t^){l-l'H^) 



K{1') - F{aTcsmT;l') 



(14) 



Together, (13, 14) invert ([9| 10) and, together with (|6| [7j), give the value of Baxter's 
parameters for equivalent points of the phase diagram. It is important to notice that 
these identities ensure that the domain (11) maps into |r| < 1 and A < — 1 and vice- 
versa. 



< 1 and Jy > 1. 



The idea of the analytical continuation is to use (|13|) to extend the domain of / 
md its 
[T| with \Jz 

Jy{uJ) 

and thus 



beyond its natural regime. To this end, let us start from a given region, let say la in fig. 

By ([T]) we have: 

cn(iti; /) 



-A 



dn(iu; I) 



JziuJ) 



-r 



dn(iu; /) 



(15) 



J? 



dt 



(16) 



(17) 



Notice that within region we can recast the second identity as 

u = K{1') + F(arcsin J^; /') , 

which makes z explicitly real and < m < 2K{V). 

We now take (16) as the definition of I and use it to extend it over the whole phase 
diagram. For instance, in region IIq (0 < < 1, < J^), using (16) we find / > 1. 
To show that this analytical continuation gives the correct results, we write / = 1//, so 
that < Z < 1 and use (A. 32) in (15), which gives 

Jy{U, I) 



dn(iM; /) 
cn(i£t; /) 



Jz{u,l) 



cn(iu; /) 



(18) 

Note that we also 



which reproduces Baxter's definitions in (6): F = —4^, A = — 

rescale the argument u = lu, so that < m < 2K{1'). This shows that we can use (15) 
to cover both regions Iq and by letting < / < oo (note that / = 1 corresponds to 
the boundary between the two regions). 

We can proceed similarly for the rest of the phase diagram. However, the second 
part of (16) needs some adjustment. Since elliptic integrals are multivalued functions 
of their parameters, to determine the proper branch of the integrand one starts in the 
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principal region and follows the analytical continuation. The result of this procedure 
can be summarized as 



lU 



K{1') + F(arcsin J,; /') + 1 - sign( Jj,) K{1) 



F(arcsin Jy] I) - K{1) + i 1 + sign( J^) K{1' 



(19) 



where all the elliptic integrals here are taken at their principal value. We did not find 
a single parametrization without discontinuities in the whole phase diagram. Close to 
~ 0, the first line of (19) is continuous, while the second has a jump. The opposite 



happens for Jy ~ 0, but, due to the periodicity properties of (|9| 10), both expressions 



are proper inversions of (15) valid everywhere. 



Thus, we accomplished to invert (15) and to assign a pair of {u,l) to each point 
of the phase diagram {Jy, Jz), modulo the periodicity in u space. The analysis of these 



mappings shows an interesting structure. From (16) it follows that that regions I's have 
< / < 1, while II's have / > 1, which can be written as / = l/l, with < / < 1. 
Finally, regions Ill's have purely imaginary I = il/l', with < / < 1. In each region, 
the argument runs along one of the sides of a rectangle in the complex plane of sides 
2K{1) and i2K{l'). Thus, we can write 

cn(^ u{z); I) , . 1 



Jyiz,l) 



JziuJ) 



< Z < TT 



where ( 



dn{C u{z); /) ' ' ' dn{C u{z); I) ' 

1,/,/' in regions I's, II's, and Ill's respectively. Moreover, u{z) = i(7r 



(20) 



\K{r) for regions of type a; u{z) = {7T-z)f^K{l) for type 6's; u{z) = 2K{l) + izlK{V 



for regions of type c; and u{z) = i2K{l') + z'^Kil) for regions of type d. In figure [2j we 
draw these four paths for the argument and the directions of the primary periods of the 
elliptic functions for regions of type I, II, and III. Figure |3] shows the contour lines of 



(u, /) in {Jy, Jz) plane, as given by (20). 

It is known that only purely real or purely imaginary arguments (modulo a half 
periodicity) ensure the reality of an elliptic function (and thus that the coupling constant 
are real and the Hamiltonian hermitian). We will show that this regular structure is 
actually connected with the symmetries of the model. The validity of this analytical 
continuation is checked, like we just did for region IIq, using the relations connecting 
elliptic functions of different elliptic parameters (see Appendix A or L21J)- These 
relations are rooted in the modular invariance of the torus on which elliptic functions are 
defined. In the next section, we are going to show explicitly the action of the modular 
group, in covering the phase diagram of the XYZ model. 

Before we proceed, we remark that, while it is possible that the validity of Baxter 
parametrization outside of its natural domain has been observed before by other author's 
analysis, we did not find any reference to it in the literature and we believe that we are 
the first to show in details how to construct the complete extension to the whole phase 
diagram, which is condensed in Table [T] in the next section. 
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• 1 
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d 











-2K{1) 

-\2K{1') 


g 2if(i) «( 

• 



Regions I 

o<;<i 



K{1) ^{K) 



Regions II 
I = l/l 
< Z < 1 



Regions III 
l = U/l' 
< [ < 1 




K{1) ^{K) 



Figure 2. On the top left, in regions of type a, b, c, and d, the argument of tlie elliptic 
functions runs along one of the sides of the rectangle. In the remaining quadrants, the 
directions of the quater-periods K{1) and iK{l') in I, II, and III regions. 




Figure 3. Contour plots of the {u, I) lines in the different regions of the phase diagram, 
as given by ( 20 1 . 
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5. The action of the modular group 



In the previous section, we showed that Baxter's parametrization of the parameters of 
the XYZ chain can be analytically extended outside of the principal regime to cover the 
whole phase-diagram and thus that any path in the real two-dimensional {Jy, Jz)-sp8ice 
corresponds to a continuous path of {u, I) in C^. This extension is due to the fact that 
the analytical continuation of elliptic functions in the complex plane can be related to 
the action of the modular group. We refer to |21) for a detail explanation of the relation 



between elliptic functions and the modular group and to Appendix A for a collection of 
useful identities. 

Points related by a modular transformations correspond to the same (F, A) point 
in the mapping to the eight- vertex model (6][7). We now want to show that the modular 
group can be connected to the symmetries of the model that we discussed in the 
introduction. Let us start again in region la, but to make the modular structure more 
apparent, this time let us write the parametrization in terms of theta functions using 



(A.6) 



03{O\r) 02 [jz - 7r)r|r] 
^2(0|r) 6, [{z - n)r\r] 



J. 



vr r r 



ir T T 



(21) 



where < ^ < tt, and we used (A.4). 



As the Jacobi elliptic functions are doubly-periodic, the torus on which they are 
defined can be viewed as the quotient C/Z^. In region la, for < / < 1, the fundamental 
domain is the rectangle with sides on the real and imaginary axis, of length respectively 
Ui = 4:K{1) and = i4K{l'), and the path followed by z is a half-period along the 
imaginary axis. A modular transformation turns the rectangle into a parallelogram 
constructed on the same two-dimensional lattice, but leaves the path of the argument 
untouched. To this end, we extend the representation of the modular group to keep 
track of the original periods and introduce the two vectors and aj, which points 
along the original real and imaginary axis. In la, these two vectors are o"/j = ^ = 1 and 
0"/ = ^ = r (we remind that, when using theta functions, every length is normalized 
by the periodicity on the real axis, i.e. uji). We will also need to extend the modular 
group and thus we introduce the additional complex number = — tto"/, which is a sort 
of an initial phase in the path of the argument. 



With these definitions, we write (21) as 
^3(0|r) [zaj 



J II 



J. 



The two generators of the modular group act on this extended representation as 



(22) 



/ . \ 



/r + l\ 



/ . \ 



<7r 



S 



(23) 



V / V 



V ^ y V y 
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A general modular transformation which gives r' 
1 



d — br' . 



also does 

a+b T 



-c + ar 



(24) 



a + hr a + br a + br 

i.e. it expresses the original real and imaginary periods in terms of the two new periods. 
For instance, applying the S transformation Sr = = — ^ to region la we have 

^3(0|r5) 02 [z - 7i\ts] ^3(0|r) 6, [{z - 7r)r|r] 



' 02iO\Ts)9s[z-n\Ts] e,iO\T)es[iz-n)T\T] ^' ^ ' 

. J ^ 03i0\rs)e, [z-n\rs] ^ 9,i0\r) 9,[iz-n)r\r] ^ 

' OMrs)Os[z-n\Ts] 9,i0\r) 9s[iz - n)T\T] ' ^^^^^ 
which is equivalent to the Baxter's elliptic parametrization of region I;, and shows that 



R 



y,z 



The other generator of the modular group Tr = = r + 1 gives 
^3(0|rr) 92 [{z - 7r)r|rT] _ ^4(0|r) 92 [{z - 7r)r|r] 



^2(O|rr)03 [{z-7i)r\rT] 
^3(01x^)^4 [{z-n)T\TT] 



92{0\t) 04 [{z-n)r\r] 
9,{0\t) 9s [iz-n)r\r] 
04(O|rr) 9s [{z - 7t)t\tt] 9s{0\t) 9^ [{z - 7r)r|r] 
which maps la into region IIIq, in agreement with ([T]) and gives 

T Ra;,2: Riy 5 

as can be seen by rescaling the couplings J^, Jy, and by 1/Jz- 
Finally, we introduce an additional operator P: 

/ . \ 



Jz 

1 

X 



(26) 

(27a) 
(27b) 

(28) 



/ r \ 



V y 



y + vr [(Tjj + 0-/] j 



(29) 



which shifts by half a period in both directions the origin of the path of the argument 
Applying it to gives 

9Mr) 92 [zr + ^\T] 



02(O|r) 9s [zt + t,\t] 
9s{Q\t) 9, [zT + n\r] 



y ' 



04(O|r) ^^3 [^r + 7r|r] 
which covers L and leads to the identification 



The three operation T, S and P, with: 



fS-T) 



I, 



I, 



I, 



P S = S P 



(30a) 
(30b) 

(31) 

P ■ T = T ■ P . (32) 
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1 
/ 
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1 
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r 

T + 1 
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r 
r+1 


1 
I 




STS 
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Table 1. List of the action of modular transformations, according to (23 29 1, in 



mapping the phase diagram, starting from with the parametrization given by (22) 



generate a Z2-extension that can be identified with PGL{2,1j). These are the same 
group laws satisfied by the generators of the symmetries of the XYZ chain, as discussed 
in the introduction. 

We collect in table [T] the action of the different transformations generated by T, 
S and P in mapping region la in the rest of the phase diagram. It is clear that & is 
isomorphic to the group generated by Rj/,2, Rx,z, and Py,z, as we anticipated, and thus, 
for instance, S • T = ^ ■ R^, ^ = Ilx,z ■ ^y,z- A few additional comments are in order. 
In the previous section we gave a prescription for the path of the argument around a 
rectangle, see (20) and figure |3} This structure is explained by table [l] the S generator 
interchanges the role of the axes and thus the path that runs along the imaginary axis 
in the new basis runs along the real one, consistently with the prescription given for 
regions of type b and d. The action of P is to shift the path and thus interchanges 
regions of type a with c and b with d (and viceversa). The role of P can be thought of 
as that of promoting the modular group to its discrete affine. 

In conclusions, we see that the whole line P G (—00,00) can be used in (20), 
which divides in three parts, one for each of the regions of type /, //, or ///. It is 
known, that, for each point on this line, only if the argument of the elliptic functions 
is purely imaginary or purely real (modulo a half periodicity) we are guaranteed that 
the parameters of the model are real. Hence the two (four) paths, related by S (and P) 
duality. The paths of the argument are compact, due to periodicity. Thus, the mapping 
between [u, I) and {Jy, Jz) is topologically equivalent to — )■ ]R(8)C(S)Z2®Z2. Extending 
the mapping by considering a generic complex argument m, would yield a non-hermitian 
Hamiltonian, but would preserve the structure that makes it integrable. We do not 
have a physical interpretation for such a non-unitary extension, but its mathematical 
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properties indicates that it might be worth exploring this possibihty. 
6. Conclusions and outlooks 

The possibility to write an exact solution for the XYZ spin chain is based on its relation 
with the eight- vertex model. However, the relation between the parameters (F, A) of 
the latter is not one to one with the couplings {Jx, Jy, Jz) of the former (and with the 
Boltzmann weights of the classical model as well): indeed such mapping changes in the 
different regions of the phase diagram, see eq.ns g jr)). In this paper, we have shown 
that this re-arrangement procedure can be reproduced by means of the natural analytical 
extension of the solution valid in a given initial region. We have also given a prescription 
that relates the parameter [u, I) of the elliptic functions and the physical parameter of 
the XYZ model in the whole of the phase diagram and provided an inversion formula 
valid everywhere. 

Clearly symmetries in the space of parameters become manifest as invariance 
properties of observables and correlations. For example, in [9] we have analytically 
calculated entanglement entropies for the XYZ model and in particular obtained an 
expression for the corrections to the leading terms of Renyi entropies that emerge as soon 
as one moves into the massive region, starting from the critical (conformal) line. Also, 
by making explicit use of such a kind of parametrization, we have shown that finite size 
and finite mass effects give rise to different contributions (with different exponents), thus 
violating simple scaling arguments. This topic deserves further analytical and numerical 
studies that we plan to present in a future paper. 

The main novelty of the approach presented in this paper lies in the connection 
between analytic continuation in the parameter space and the action of the modular 
group. This connection allowed us to show that a certain extension of the modular group 
realizes in parameter space the physical symmetries of the model. This symmetries are 
in the form of dualities, since they connect points with different coupling, related by 
some spin rotation. 

Modular invariance plays a central role in the structure and integrability of 
Conformal Field Theories in 1 + 1 dimension. On its critical lines, the XYZ model 
is described by a c = 1 CFT and thus, in the scaling limit, its partition function is 
modular invariant in real space. Our results show that, even in the gapped phase, it 
is also a modular invariant in parameter space, due to the symmetries of the model. 
Elliptic structures are common in the solution of integrable models and it is tempting 
to speculate that their modular properties do encode in general their symmetries and 
thus the class of integrability-preserving relevant perturbations that drive the system 
away from criticality. 
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Appendix A. Elliptic Functions 

Elliptic functions are the extensions of the trigonometric functions to work with doubly 
periodic expressions, i.e. such that 

f{z + 2nuJi + 2mw3) = f {z + 2nuji) = f{z) , (A.l) 

where 0Ji,Us are two different complex numbers and n,m are integers (following |21] . 
we have ui + ijJ2 + cos = 0) . 

The Jacobi elliptic functions are usually defined through the pseudo-periodic theta 
functions: 



9,{z- q)=2 ^(-l)"g("+i/2)2 ^-^^^^^ ^ ^ 2a) 

n=0 

oo 

92{z; g) = 2 ^ qin+m' pQg[(2n + l)z] , (A.2b) 

n=0 

oo 

^3(2;g) = 1 + 2^g"'cos(2n2) , (A.2c) 

n=l 

oo 

0^{z; q) = l + 2 cos(2n2) , (A.2d) 

n=l 

where the elliptic nome q is usually written as q = e™"^ in terms of the elliptic parameter 
T and accordingly 

e^{z\T)^e^{z-q). (A.3) 
The periodicity properties of the theta functions are 

6li(2|r) =- 6*1(2 + 7r|r) = -A6'i(2 + 7rr|r) = A6'i(2 + vr + 7rr|r) , (A.4a) 

6l2(2;|r) =- 6*2(2; + 7r|r) = A6'2(z + 7rr|r) = -A6'2(2; + vr + 7rr|r) , (A.4b) 

Oz{z\t)= e^{z + 'K\T)= \e^{z + -KtIt) = \e?Xz + -K + -KtIt) , (A.4c) 

94{z\t) = 614(2; + 7r|r) = -A6'4(z + 7rr|r) = -A6'4(z + TT + 7rr|r) , (A.4d) 
where A = ge^'^. 

Moreover, incrementation of z by the half periods ^tt, ^ttt, and |7r(l + r) leads to 

9i{z\t) = — 92{z + ^7r|r) = —i^9i{z + |7rr|r) = —i^9z{z + |vr + |vrr|r) , (A. 5a) 

92{z\t)= 9i{z+\'k\t)= h9^{z + \'kt\t) = ^Ji,9i{z+\TX + \'kt\t) , (A.5b) 

9^{z\t)= 9i{z+\Ti\T) = ^92{z + \'kt\t) = ^i9i{z + \ti + \'kt\t) , (A.5c) 

9^{z\t)= ^3(2;+ ivrlr) = -i/i^i(2;+ |7rr|r) = i/x^2(^ + + iyrrlr) , (A.5d) 
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with fi = g^/^e^^. 

The Jacobi elhptic functions are then defined as 



sn(u; k) 
cn{u] k) 
dn('u; k) 



where u = 91{{]\t) z 
K{k) = 



^3(0|r) 


0i{z 


r) 


^2(0|r) 




T) 


^4(0|r) 


e2{z 


t) 


^2(0|r) 


e,iz 


T) 


OMr) 


e-siz 


r) 


^3(0|r) 


e^iz 


T) 


lK{k)z, 


T = 


.K(k') 



k' 



dt 



a/1 — k^, and 
d^ 



1 - A;2 sin^ 6 



(A.6a) 
(A.6b) 
(A.6c) 

(A.7) 



/o ^{l-t^){l-kH^) 

is the elhptic integral of the first type. Additional elliptic functions can be generated 
with the following two rules: denote with p, q the letters s,n,c, d, then 

1 , , s pn{u; k) 



pq(M; k) = 



pq{u; k) 



(A.8) 



qp('u; k) ' ^^'^ ' ' qn(M; k) 

For < < 1, K{k) is real, iK'{k) = iK{k') is purely imaginary and together they 
are called quater-periods. The Jacobi elliptic functions inherit their periodic properties 
from the theta functions, thus: 

sn(M; k) = - sn(M + 2K; k) = sn{u + 2iK'; k) = -sn{u + 2K + 2iK'; k) , (A.9a) 
cn(M; k) = - cniu + 2K; k) = -cn{u + 2iK'; k) = cn{u + 2K + 2iK'; k) , (A.9b) 
dn(M; k) = dn{u + 2K; k) = -dn(u + 2iK'; k) = -dn{u + 2K + 2iK'; k) . (A.9c) 



The inverse of an elliptic function is an elliptic integral. The incomplete elliptic 
integral of the first type is usually written as 

F((f);k)= = = sn-^(sin0;fc) , (A.IO) 
' ^l-Psin^^ 



and it is the inverse of the elliptic sn. Clearly, = K{k). Further inversion 
formulae for the elliptic functions can be found in [21j. 
Additional important identities include 

k - k(r) - ^ k' - k'(r) - ^ (A 11) 

^ ~ ~ ' k -k[r)- ^^,^^1^^ , (A.ll) 

and 



sn^{u]k) + cn^{u]k) =1, (A.12) 
dn^{u] k) + k'^sn^{u] k) = 1 , (A.13) 
dn2(M; k) - k\n\u] k) = k'^ . (A.14) 

Taken together, the elliptic functions have common periods 2uJi = 4K{k) and 200^ = 
4:iK'{k). These periods draw a lattice in C, which has the topology of a torus. While the 
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natural domain is given by the rectangle with corners (0, 0), (4i^, 0), (4/^, 41/^'), (0, 4iJC'), 
any other choice would be exactly equivalent. Thus, the parallelogram defined by the 
half-periods 

uj'-^ — aui\-\-huiz-i ui'^^ — cuii + duizi (A. 15) 

with a, 6, c, d integers such that ad — he — 1, can suit as a fundamental domain. This 
transformation, which also changes the elliptic parameter as 

, cjo c + dr 



' ~ ~ ~ — TIT ' 
can be casted in a matrix form as 



I 



v 



a b 
c d 



\ ( \ 



(A.16) 



(A.17) 



and defines a modular transformation. The modular transformations generate the 
modular group PSL(2,Z). Each element of this group can be represented (in a not 
unique way) by a combination of the two transformations 

/ 



S = 



T = 



1 

^ 0^ 



T = 



1 
T 



V 



1 1 



r' = T + 1 , 



which satisfy the defining relations 



{STf 



The transformation properties for the S transformation are 

1~ 



02 ( z\ 



e. z 



— i(iT) 2e 6i{tz\t) , 
(— ir)2e Q/^\tz\t) , 

1 \TZ^ 

(— ir)2e ^ Qj,[tz\t) , 



1 \TZ^ 

-it) 2 6^^6*2 fr^; I r) 



Prom which follows 



- A;'(r) , 



A;' 



(A.18) 
(A.19) 

(A.20) 

(A.21a) 
(A.21b) 
(A.21c) 
(A.21d) 

(A.22) 



and 



X l^-i^ = K'{t) = X(A;') , \K' (^-^^ = \K{t) = \K{k) . (A.23) 
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, sn(m; k) 
cn(m; k) ' 
dn(m; k) 



cn(M; k') 



cn(m; /c) 



cn(m; A;) ' 
For the T transformation we have 

^i(^|t + 1) = e'^'%{z\T) 

e^{z\T + i) ^e^{z\T) , 



and thus 



and 



Finally 



K{t + 1) = k'{T)K{T) , 

sn{u/k'; k) 



k'{T + 1) 



A;'(t) 



(A.24) 

(A.25a) 
(A.25b) 
(A.25c) 
(A.25d) 

(A.26) 



sn m; 1 



k 



iK'{T + l) = k' [K{t) + iK'ir)] . (A.27) 
cn{u/k'; k) 



cn I 1 



dn{u/k'; k) 



(A.28) 



dn(M//c'; k) ' 
^dn(VA;';A;) ' 

While it is true that by composing these two transformations we can generate the 
whole group, for convenience we collect here the formulae for another transformation 
we use in the body of the paper: r — > j^^, corresponding to STS 

iiF^i((l-r)^|r) , 
Fe^({l-T)z\T) , 
F^2((l-r) 



02 (z\ 



04 ( Z 



r 

1 -T 

r 

1 -T 

r 



1 -r 

r 

1 - r 



u r 



^ =i^F^4((l-T)z|T) , 



(A.29a) 
(A.29b) 
(A.29c) 
(A.29d) 



where F = (1 — T)2e'' . We have 



r 



1 -T 



k{T) k ' 



A;' 



r 



1 -T 



i^=i^, (A.30) 



k{r) k 



and 
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Moreover, 



sn ( m; - I = k sn ( k] , cn ( u: - I = dn ( k 



k J V/c'/' \ ^ k j V/c' 



dniu;]-] =cn(j-k] , ... (A.32) 
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